Nonlinear conductivity of diffusive normal-metal contacts 
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Metal microbridges with a high impurity content and shorter than the energy relaxation length 
are considered. Their conductance is calculated with allowance made for the Coulomb electron- 
electron interaction. It is shown that nonequilibrium electrons in the microbridges gives rise to a 
nonlinear current-voltage characteristic. 



It is common knowledge that tunnel junctions between 
disordered metals show nonlinear conductance at low 
temperatures (the so-called zero-bias anomalW), which 
results from the electron-electron interaction.EJ On the 
other hand, short microbridges possess meny properties 
of tunnel junctions. Therefore, it is interesting to inves- 
tigate their nonlinear properties at low temperatures. 

Nonlinear conductivity of weakly disordered metals at 
lew., temperatures was, investigated both experimental- 
lyHtl and theoreticallyu in a number of papers. In these 
papers, the nonlinear behavior resulted from the heat- 
ing of the electron gas by the current and from the 
temperature-dependent correction to the conductivity 
arising either from weak localization or electron-electron 
interaction effects. The resulting current-voltage char- 
acteristics were influenced, however, by several param- 
eters describing the heat transfer from the electron gas 
in a massive sample, e.g., by the electron-phonon cou- 
pling constant and the acoustic transparency of the sam- 
ple boundary. Therefore it was difficult to compare the 
theoretical results with theoretical predictions. For metal 
microbridges, the situation is different because the heat 
transfer is determined by the diffusion of hot electrons 
into the banks and no additional parameters are required. 

In this paper we investigate the influence of Coulomb 
exchange electron-electron interactions on the current- 
voltage characteristics of disordered-metal microbridges 
shorter than the energy relaxation length. The char- 
acteristic energy scale in this case is the largest of the 
quantities eV and kT, wher V is the bias voltage and 
T is the temperature. Therefore the electron distribu- 
tion function is essentially nonequilibrium at low tem- 
peratures kT <C eV, and one may expect a nonlinear 
current-voltage dependence. 

Consider a metal microbridgc connecting two massive 
banks with the x axis directed along the contact. The mi- 
crobridge length L is assumed to be larger than the elec- 
tron elastic mean free path and the microbridge width. 
Yet, L is assumed to be shorter than the energy relax- 
ation length. 

In studying the current-voltage characteristics of the 
contact under nonequilibrium conditions, the Keldysh di- 
agrammatic techniqueO may be conveniently used. In this 



technique, three types of Green functions are used, 
G R (1, 2) = i6(t 2 - t 1 )(^{\)^+{2) + ^+(2)^(1)), 

G A (l, 2) = -<0(ti - t 2 )(V(l)V + (2) + V+WCl)}, 

G K (1, 2) = -i(V(l)V + (2) - ^+(2)^(1)). 

In the absence of the electron-electron interaction, the 
impurity-averaged retarded and advanced Green func- 
tions are given by the well-known expressions 



G R ( A \e,r, br') 



d 3 p exp[ip(r — r')] 
(2tt) 3 e -p 2 /2m±i/2T' 



(1) 



where r is the elastic relaxation time. The Green func- 
tion G K obeys the integral equation 

G*(£,r,r') = ^^ J d 3 n G«( £ ,r,n) 

xG K {e,vx,T l )G A {e,v i y), (2) 

where TV is the electron density of states at the Fermi 
level. By setting r = r' in (||), one may obtain a dif- 
fusion equation for the Green function G K (e, r,r') (see 
Ref.EI) 



£>V 2 G A >,r,r) = 0, 



D = -vir. 



(3) 



In the case of massive banks, the electron distribution at 
the ends of the microbridge is equilibrium and therefore 
the boundary conditions for Eq. (S) are 



G K (e, r, r)\ x= _ L/2 = -2ttN[1 - 2n(e - eV/2)}, 



G K {e,v,v)\ 



x=L/2 



-2tvN[1 - 2n(e + eV/2)}, (4) 



where n{e) is the Fermi distribution function and V is 
the voltage drop across the contact. The current flowing 
through the contact may be expressed in terms of the 
Green function G K via the formula 

3w = -— lde(^--4i\ G k {s,t,t')\ ,. (5) 
J 2m \dx dx' K ' lr=r ' 



1 




FIG. 1. The diagrams describing the current to first order 
of the electron-electron interaction. 

The corrections to the current from the electron-electron 
interaction are represented by the three diagrams shown 
in Fig. 1 and their complex conjugates. In these di- 
agrams, solid lines denote the advanced and retarded 
electron Green functions, the black circles denote the 
Green functions G (r, r). The single dashed lines repre- 
sent the impurity-potential correlator (27riW) _1 5(r — r'). 
The shaded rectangles denote the impurity-averaged two- 
particle Green functions 

P(u, r, r') = (G A (e + u, r, r')G R (e, r', r))„„ p , (6) 

which obey the equation 

(iuj + DV 2 )P(uj, r, r') = -2nN5(r - r'), 

P( W ,r,r')L =±L/2 =0. (7) 

The shaded semicircles denote the impurity-renormalized 
electron vertex (2ttNt)~ 1 P(oj, r, r'). The electron- 
electron interaction is represented by the wavy lines. 
Taking into account the Debye screening, one obtains 
the following equation for the retarded potential of the 
electron-electron interaction, 

DV 2 V r (uj, r, r') = JV _1 (-ta; + DV 2 )5(r - r'), 

V R ( W ')\ x=±L/2 = 0. (8) 

The advanced potential V A is the complex conjugate of 

yR 

These diagrams result in a contribution to the current 
density in the form 

ji(r) = 2(2n)- 5 eDN- 3 

xlm J de J duj J d 3 r x J d 3 r 2 G K (e,r,r) 
xP(-lo, r, v 1 )G K {e - w, r U Ti)V R (w, r u r 2 ) 

xV r P(- W ,r 2 ,r). (9) 

To take into account the current-conservation law, this 
current density should be averaged over the contact 
length. 

Substitute the solutions of (0), and (||) into (|). 
In the low-temperature limit kT <C eV, the Fermi 
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FIG. 2. The differential conductance versus voltage. The 
inset shows the electron distribution function at an arbitrary 
point inside the contact. 

distribution function is step-like and the integration over 
e may be easily performed. The integration with respect 
to u> in (|^) may be conveniently replaced by the integra- 
tion with respect to a dimensionless quantity ?/, which is 
related to u> by y 2 = (L 2 /2D)lo. Therefore, the correc- 
tion to the current flowing through the contact may be 
presented in the form 

where s = L(eV/2hD) 1 / 2 and 

Q(y) = {(l/4Msmh(2y) - sm(2y)] 

-y _1 [sinh(2?/) + sin(2y)] 

+2y~ 1 [smh(y) cos(y) + sin(y) cosh(j/)]} 

x [cosh(2y) -cos(2y)]" 1 . (11) 

The corresponding differential resistivity versus voltage 
curve is shown in Fig. 2. At high voltages eV 3> TiD/L 2 , 
h w -(8/n)(DeV/2hL 2 ) 1 / 2 . Therefore, the voltage 
dependence of resistivity in this voltage range has the 
shape characteristic of the temperature dependence of 
resistivity in one-dimensional conductors.EJ This is quite 
natural, because the width of the drop in the electron 
distribution function is now determined by the applied 
voltage and not by the temperature. Therefore, the 
characteristic momentum transfer in (^|) is of the order 
of (eV/hD) 1 / 2 . In the opposite case of small voltages 
eV <C UD I L 2 , the momentum transfer is limited by %/L. 
Therefore, the conductivity tends to a constant value, 
i.e., 1\ = — 0. 17 (e 2 /Ti)V. Hence, increasing voltage gives 
rise to a cross-over from a zero-dimensional case to a 
quasi one-dimensional case. 
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Note that the effect of the electric field is not reduced 
to a simple heating of the electron gas to some effective 
temperature. Instead, the electron distribution function 
consists of two subsequent steps positioned at ef ± eV/2 
(see the inset in Fig. 2) and therefore has an essentially 
non-Fermian shape. In principle, the noncquilibrium 
electron distribution may also affect the phase-breaking 
time and therefore the weak localization contribution to 
the conductances However, weak localization corrections 
may be suppressed by a sufficiently strong magnetic field, 
and the universal conduction fluctuations&Q may be el- 
eminated by averaging the results over a number of sam- 
ples. 
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* Very recently, the nonlinear Aronov - Altshuler correction 
to the conductance of diffusive mesoscopic contacts was ex- 
perimentally studied by H.B. Weber et al., Phys. Rev. B 63, 
165426 (2001). Though not cited there, these theoretical re- 
sults are relevant to this problem. 
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